RPA GREEN'S FUNCTIONS OF THE ANISOTROPIC HEISENBERG MODEL 

Andre Johannes Stoffel and Miklos Gulacsi 
Max-Planck-Institute for the Physics of Complex Systems, D-01187 Dresden, Germany 
Nonlinear Physics Centre, Australian National University, Canberra, ACT 0200, Australia 

November 10, 2009 

We solve in random-phase approximation the anisotropic Heisenberg model, including 
nearest and next-nearest neighbour interactions by calculating all Green's functions and pair 
correlation functions in a cumulant decoupling scheme. The general exposition is pedagogic 
in tone and is intended to be accessible to any graduate student or physicist who is not an 
expert in the field. 



I. INTRODUCTION 



Recently we analysed the properties of an anisotropic Heisenberg model in an external 
longitudinal field on a bcc lattice, with a particular application to supersolids 
However, detailed derivation of the Green's functions and correlation functions have not yet 
been published, as such we fill this gap hereafter. For the interested reader we present in 
Appendix A the connection between the model analysed hereafter and 4 He. 

We work with the standard anisotropic Heisenberg model, defined by the Hamiltonian: 

i ij ij 

on a bcc lattice shown in Figure ([IT) . It is known that the Hamiltonian defined in Equation ([1]) 
gives four magnetic phases: ferromagnetic, canted ferromagnetic, canted anti-ferromagnetic 
and anti- ferromagnetic phases which we abbreviate by FE,CFE, CAF and AF, see Table (JI]). 
The order parameters, mi for off-diagonal long range order and m<i for diagonal long range 
order in the magnetic system are defined by: 

mi = (S% + (£%) 

m 2 = (S Z A ) - (S Z B ) (2) 



In the following we will use these order parameters to identify the phases within the random- 
phase approximation. 
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FIG. 1: The bcc lattice consists of two interpenetrating sc sub-lattices, i.e., sub-lattice A and 
sub-lattice B. For simplicity we only have drawn the two dimensional case. 



The Heisenberg model has been studied thoroughly although there exists only a classical 
mean-field solution for the canted ferromagnetic and anti- ferromagnetic phases 5|,|6|]- Hence, 
our interest is in obtaining a better than mean-field solution of these canted phases. We 

nnri 

used [2|, |3j4j the equation of motion technique for the time-temperature-dependent Green's 
function 7|, G^J (t)R et /Adv = Tid(±t)([Sf(t), Sj})}, for all four phases, which we solve in 
random-phase approximation (RPA). In order to maintain pair correlations as accurate 
as possible, we have chosen a cumulant decoupling scheme to approximate higher order 
Green's functions beyond the random-phase approximation. This is the first time that a 
thorough cumulant random-phase approximation many-body calculation has been solved for 
the canted phases and as such it represents a vital step in the understanding the behaviour of 
the three dimensional anisotropic Heisenberg model in the presence of an external magnetic 
field. Details of the Green's functions are given in the next Chapter. 

By using these Green's functions we calculate exactly all the pair correlations, also given in 
the next Chapter, where we show that this decoupling scheme couples six Green's functions, 
one for each spin component in x, y and z direction on the two sub-lattices respectively, 
to a set of six equations. Due to the enormous number of terms (1024 in total) within 
these Green's functions we have decided to reproduce their exact form only at the end. 



Using these Green's functions it can be shown or alternative 
theorem of gapless modes imposes an additional condition 



argued that the Goldstone 
on the mean fields of the 
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Spin Configuration 


Magnetic Phase 


ODLRO 


DLRO 


TT 


FE 


No 


No 


// 


CFE 


Yes 


No 


// 


CAF 


Yes 


Yes 


Tl 


AF 


No 


Yes 



TABLE I: Possible magnetic phases of the anisotropic Heisenberg model. All phases are defined by 
their long range order, i.e., off-diagonal long-range order (ODLRO) and diagonal long-range order 
(DLRO). 

canted phases, reducing the number of order-parameters by two. As this condition does 
not apply to the ferromagnetic and anti-ferromagnetic phases, their Green's functions are 
structurally different. 

For completeness we also re-derive the classical mean-field Green's functions in Appendix 
B. While in Appendix C we establish a link between the classical mean-field and the random- 
phase approximations. 



II. GREEN'S FUNCTIONS 



The most versatile and successful method of solving_many-body problems involve the 
concept of Green's function. The Matsubara formalism 7| for imaginary time Green's func- 
tions at finite temperature is commonly used to analyze fermionic and bosonic many-body 
systems. However, in spin systems a different type of Green's function is primarily applied 
as this makes it easier to deal with the algebraic properties of spin operators. The retarded 
and advanced Tyablikov [9|, H0[ commutator Green's function at finite temperature defined 
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in real time are: 



G f SLW=<fH)<[5 r r(*) > ^]> (3) 

The average () involves the usual quantum mechanical as well as thermal averages and /i 
and v are elements of {x, y, z} and % and j denote the lattices sites. There is no Wick's 
theorem for spin systems and therefore no perturbative approaches by means of Feynman 
diagrams are available. However this is not necessarily a disadvantage as the method of the 
equation of motion is equally powerful in many cases and often much more compact. The 
basic idea of this method is to find a linear differential equation for the Green's function. 
Therefore we can differentiate the xy-Green's function with respect to time: 

idtGiLM = S(t)([Sf,S^)-i9(t)([[S!,H],S^) 
id t G iLJt) = 5(t)([S?,SV])+ie(-t)([[S?,H],S?\) 

(4) 

Here we have used the Heisenberg equation of motion id t S x = [Sf,H] for the operator S x . 
For the anisotropic Heisenberg model the commutator yields: 

[Sf, H] = -inSf -2iJ2 4 S i S ! + 2% E J U S - S j ( 5 ) 

3 3 

Consequently the right hand side of Equation (j4j) involves higher, third order Green's func- 
tions. In principle those third order functions can be expressed in terms of even higher order. 
This procedure will go on indefinitely yielding a series of infinite order. The central idea is 
to introduce a suitable approximation which truncates this series. Here, in order to keep the 
complexity to a manageable level we split up the three operator correlation functions. 



The cumulant decoupling 11] is based on the assumption that the last term of the 
following equality is negligible: 

(ABC) = 

(A)(BC) + (B)(AC) 

+ (C)(AB)-2(A)(B)(C) 

+((A-(A))(B-(B))(C-(C))) (6) 
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This approximation is justified if the quantum fluctuations are small and do not deviate far 
from their mean field values, which is the case in three dimensions [12]. As a consequence 
the third order correlation functions split into product terms of a second order correlation 
functions and the mean-field of a single operators: 

m(t),s^t)},s]])^ 

(Sm([S P k(t),S»\) + (S p k (t))([S?(t),S]\) (7) 
Finally we obtain a differential equation which only involves second order Green's functions. 

id t G^{t) = 

i 

+2* Jlm{t))G^{t) + (Sf(t))Gt](t)) (8) 
i 

Here we have dropped the subscripts for the retarded and the advances Green's functions 
since the equations equally hold for both. In the next Chapter we will see that the advanced 
and the retarded Green's functions are actually represented by the same Fourier transform; 
the difference comes from the path along which the inverse Fourier integral is carried out. 
Similar relations also exist for the time derivatives of G^J(t) and G^J(t). Together these 
equations form a closed set of linear differential equations. Usually linear differential equa- 
tions are readily solved by Fourier Transformation. The time-frequency Fourier and inverse 
Fourier transforms are given by: 

= J due**G°g(u>) (9) 
<W = ±-J dte**C%(t) (10) 

One might think that due to broken translational symmetry of the canted anti-ferromagnetic 
and the anti-ferromagnetic phases, featuring diagonal long-range order, a spacial Fourier 
transform might not be applicable. Actually, the canted anti-ferromagnetic and anti- 
ferromagnetic states do exhibit discrete translational symmetry, namely they are invariant 
under — > + n\a\ + ^2(12 + ^3^3, where a\, (12 and are the basic lattice vectors of 
the sc sub-lattice and ni,ri2,n3 G Z . This translation maps each sub-lattice onto itself and 
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the number of equation doubles as we have to treat each sub-lattice separately. The Fourier 
transform into k-space is defined by: 



G%(w) = J d 3 ke ik ^G xy (k,u) 



(11) 



After successively carrying out time and space Fourier transforms we derive a set of six 
algebraic equations, determining six Green's functions that represent the xy-,yy-, and zy- 
spin correlations on each sub-lattice. The detailed calculation is carried out at the end of 
this Chapter. In matrix form this set of equations reads: 



M ■ 7 



where 



M = 



( 







Ml 3 


M14 


Mis 


M 16 \ 







iuj 


M 23 


M 24 


M 25 


M 26 




-Ml 3 


-M u 


iuj 





M 35 


M 36 




-M 23 


-M 24 





iu; 


M 45 


M 46 




M 51 


M 52 


M 53 


M 54 







V 


M 61 


M 62 


M 63 


M 64 








and 
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\ 



as well as 



(27T) 



V ^» y 

(S' B ) 



V <-ss> J 



(12) 



(13) 



(14) 



(15) 



The components of the matrix M are given by: 

M 13 = 2h z + A{S z A )(4(0) - Jj(k))+A(S Z B )4(0) 
M 14 = -A(S Z A ) Jj(k) 
M 23 = -A(S z B )jJ(k) 

M 24 = 2h z + 4(^)(J|(0) - jj(k)) + A(S Z A )4(0) 
M 15 = HS« A )(4(k) - J 2 T (0)) - 4(^)^(0) 
M 16 = A(S y A )4(k) 
M 25 = A(S y B )4(k) 

M 26 = 4<S»)(j}(*) - J 2 T (0)) - 4(^)^(0) 
M35 = -MS* A )(4(k) - J 2 T (0)) +A(S%)Jj(0) 
M 36 = -A(S x A )4(k) 
M A5 = -A(S x B )4(k) 

M m = -A(S x B )(4(k) - Jj(0))+A(S* A )JJ(0) 
M 51 = 4(^)^(0) + 4(^)(J 2 T (0) - Jj(k)) 
M 52 = -A(S y A ) Jj(k) 
M 61 = -A{S B )jJ(k) 

M 62 = A(S y A )jJ(0) + A(S B }( J 2 T (0) - Jj(k)) 
M53 = -4(^)^(0) - 4(^)(J 2 T (0) - Jj(k)) 
M u = A(S x A )jJ(k) 
M m = A(S x B )jJ(k) 

M 64 = -4(S%)J?(0) - 4(^)(J 2 T (0) - J 2 T (k)) 

(16) 

Here the k-dependent coupling constants are defined by Jj(k) = J 1 T 7i(/c), Jj(k) = Jjl2(k), 
Jj(k) = Jj^i(k) and Jj(k) = Jj^k), where 

7l (k) = -X^ e ikaAB 

72 (k) = I^e ika - (17) 

O-AA 
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are the lattice generating functions. On the bcc lattice these lattice generating functions are 
given by 

7l (k) = cos ^— j cos^—j COS j 

. cos (k x ci) cos (k v a) cos (k z a) , . 

72 (k) = — — ^+ — ( is ) 

where a is the lattice constant of a simple cubic sub-lattice. Again, the spontaneously broken 
U(l)-symmetry gives us the freedom to set (S A ) = (S B ) = 0. This reduces the number of 
non-zero matrix components of M. 

It has been shown that the commutator Green's functions must not have a zero frequency 
pole 8j. This results directly from the fact that the commutator of to spin operators at long 
time distances becomes zero: \im t ^ 00 ([S tJ '(t)S 1 '}) = 0. In the present calculation the Green's 
functions do actually acquire a zero frequency pole, as given by the eigenvalues of the matrix 
M. In the ferromagnetic and anti-ferromagnetic phases this zero frequency pole is readily 
canceled out, but for the canted ferromagnetic and the canted anti- ferromagnetic phases this 
imposes an additional constraint: 

h* + 2(S Z A )(4 - Jj) + 2{S Z B )4 = 2J 1 T p!<£l) 

h* + 2{S Z B ){JI - Jj) + 2{S Z A )J\ = 2J 1 T p|(^) (19) 

(20) 

Note that these two conditions are identical to the classical mean-field equations aaa. 
We use these two relations to replace the external magnetic field h z in the matrix M. To 
combine terms we introduce following variables: 



A 1 = j 2 '(i_ 72 (fc)) + ^ 
A 2 = Jj(l - l2 {k)) + Jl 



(SI) 
(S%) 



B 1 = jJ(l- l2 (k)) +Jl { ^§{ 

-4(4'" Jj)l2 (k)(S* A ) 2 
B 2 = jj(l- I2 (k)) + Jj { ^ 

-A(4-jJ) l2 (k)(S%) 2 
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D 



-^i(k)(Jj{S* A ){S z B ) + 24{S%){S%)) 



(21) 



Then the matrix M in the canted anti-ferromagnetic and canted ferromagnetic phases, where 
Equation (fT9|) holds, is given by (in block-form): 



Mr. 



where the 2 x 2-blocks are given by: 



( Mfj Mf 2 
M£ M 2 C 2 M 2 C 3 
M 3 C ! M 3 C 2 M 3 C 3 



M 
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m; 



22 



M: 



33 



M 



12 



-M 



21 



iu 
zo; 

4^(5-1) -4C(S' A ) 
-4C(S' B ) AA 2 (S%) 



M 



13 



M: 



31 








m: 



23 



32 



4Ai(S!) 2 -Bi 4C(5^)(5^)+D 

4C(S^)(S|,)+D 4A 2 (gj,) 2 -B 2 
(51) (S|> 

4 -A (51) 4C(5i) 
4<7<S|> -4A 2 (S|) 



(22) 



(23) 



In the phases that do not exhibit off-diagonal long range order, namely the ferromagnetic 
and anti-ferromagnetic phases, Equation ( 1191) does not hold. However, as (S x ) = (S y ) = 
the matrix M reduced to a matrix of 4 x 4-dimensions: 



Afi 3 M u \ 



iu M 23 M 24 
-M 13 -M M zcj 
\ -M 23 -M 24 y 

Equation (fl2|) readily implies that the Green's functions are given by: 



(24) 



M~ 



■ v 



(25) 
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As these Green's functions in explicit form are rather complex, we defer presenting them in 
their full form until the end of this Chapter. Those readers who are not interested in the 
details of these functions might safely skip the corresponding section of the Chapter. 

The Green's functions derived above determine the state of the system and all relevant 
macroscopic and thermodynamic properties can be extracted from them. Nevertheless the 
Green's functions as derived in this random phase approximation are function not only of 
the external field h z but also of the mean fields of the spins. Therefore we have to define 
self-consistency equation which determine those mean-fields. 

III. SELF-CONSISTENCY EQUATIONS AT ZERO TEMPERATURE 

The Green's functions as derived in the previous Chapter do not bear an explicit de- 
pendence on the temperature but rather depend on the temperature through the various 
mean-fields. Thus, upon deriving the determining self-consistency equations the tempera- 
ture will be introduced explicitly into the formalism. We have decided here to separate the 
zero temperature and the finite temperature formalism. Even though, it would generally be 
easy to carry out the limiting process T — * at any time it would cause numerical difficulties 
to do so at a later stage as the temperature usually appears in terms of 1/T. Also we find 
it instructive to derive the zero temperature formalism separately as it provides unveiled 
insight into the quantum nature of the model. 

Naturally there seem to be three ways to set up self-consistency equations. We could, for 
example calculate: 



an equivalent relation holds for (Sf) and G*j. Here, for the sake of readability we have dis- 
regarded the sub-lattice subscripts A and B. Unfortunately, these equations are not suitable 
to calculate the mean-fields self-consistently but rather lead to an identity, giving identically 



are exactly those ones contained in the right hand side of Equation (fl2|) . Therefore the 



5 (5?> = <i[fi?(0),5?]> 
= G% (O)-GS' (0) 




(26) 



\(S Z ) on both sides. That is due to the structure of the equation of motion, and the (Sf) 
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FIG. 2: The correlation functions are given by contour integrals which circumscribe the four poles. 
The residue theorem states that the integral equals the sum over the four residues. 



appropriate and only choice to define self-consistency equations is: 



(27) 



In order to establish the link between the correlation functions and the corresponding Green's 
function the following spectral expansion of the Green's function at absolute zero as is readily 
obtained from Equation (J3]): 

(n \Sy\n){n\S?\n ) 



l JRet/Adv y ' 



E 



UJ — [U)~ 



uj no ) ±ie 



E 



(n \SV\n){n\S?\no) 
uj + (uj n - uj nQ ) ± ie 



(28) 



where we have used that: 



1 roc 

2vr 



6(t)dt = lim 



(29) 



ho 27r uj — uj 1 + ie 

Here |n ) is the ground state of the system and \n) refers to the complete set of eigenstates. 
The spectral representation (Equation ( J28l) ) shows that the integral: 



/oo 
dwe^XGH (uj)-Gff (uj)) 
y l JAdv y ' l jRet y >> 

■OO 

due^(G%(uj + ie)-G%(uj-ie)) 



lim 



(30) 



is a contour integral enclosing the poles of the Fourier transformed Green's function, as 
can be seen in Figure ([2]). The path in the upper and lower half planes correspond to the 
retarded and Green's functions respectively. According to the residue theorem of complex 
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analysis the value of the integral is given by the sum of the residues of the enclosed poles. 
As u no is by definition the lowest energy state, the spectral representation of the Green's 
function shows that the correlation function (Sf(t)Sj) corresponds to all negative poles and 
and the conjugated correlations function (S y Sf(t)) corresponds to all positive poles. Hence 
we obtain the desired correlation function (Sf(t)S y ) by restricting the contour integral to 
negative frequencies: 

-i lim / e wt G yy (u + ie) - G\?{u - ie)du 

€ r ® J QO 

+(n \Sf\n ){n \S]\n ) ( 31 ) 

The second term on the right side appears because its contribution in the spectral represen- 
tation (Equation (|28j) ) is canceled out. As we have broken the U(l) symmetry of the ground 
state in a way that (S y ) = this term will yield zero anyway and we will disregard it in 
further discussion. Using the Green's function's Fourier transform into k-space we obtain: 

-i lim J J d 3 k due i{ - kr ^ +ut) [G yy (k, u + ie) 

-G yy {k,u -ie)]dw 
= 2it J d 3 ke ihr 'i [Residue (G yy (k, u)e iuJ \ -lo x ) 

+Residue (G yy (k, u)e iuJ \ -u 2 )\ 

(32) 

Now we can define two self-consistency equation which determine the spin fields: 

F A := 2tt J d 3 fc[Residue(G^(fc,a;),-c<;i) 
+Residue(G^(fc, u), -u 2 )\ - - = 



F B := 2tt J rf 3 fc [Residue(Gf (k,^),-^) 



+Residue(Gf (Jfe, u), -u 2 )} - - = (33) 



In the canted anti-ferromagnetic and canted ferromagnetic phases we use Equation ([19]) to 
eliminate (S^) and (S%) and re-introduce h z instead and obtain: 
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F c B ({S* A ),{S%),h z ) = 



(34) 



similarly, we obtain the self-consistency equations for the ferromagnetic and anti- 
ferromagnetic phases: 



the three dimensional numerical integral contained in those equations can be reduced by 
one dimension by introducing a two dimensional density of states pij'i,^) = J cPkd^^ — 
7i(fc))5(72 - 72(&)), so that: 



Consequently, in the self-consistency equations the DOS p(7i, 72) is the only remaining term 
depending on the specific lattice structure. (In the actual numerical calculation it is, in or- 
der to avoid singularities in the origin, feasible to integrate over p(arccos(7i), arccos(72)).For 
small k, arccos(7i) and arccos^)) vary as ~ k and therefore take the important contribu- 
tions at small k more accurately into account than 71 and 72 which vary as ~ k 2 .) Through 
this procedure we have archived a wider applicability of the corresponding equations as other 
systems, i.e. lattice geometries exhibiting frustration or two dimensional systems where lin- 
ear spin-waves are still a valid approximation, can easily be accomplished by simply inserting 
the appropriate DOS. The DOS for the bipartite bcc lattice is shown in Figured 

IV. SELF-CONSISTENCY EQUATIONS AT FINITE TEMPERATURE 

At finite temperature the system additionally exhibits thermal fluctuations expressed 
thorough averaging with the Boltzmann weight. As a consequence the negative and positive 
poles of the Green's functions do no longer refer solely to the corresponding correlations 
functions or their conjugates respectively. As the temperature is turned on the contribution 
of the poles to the correlation functions starts to overlap with the contribution of their con- 
jugate correlation functions, where the degree of overlap is determined by the Boltzmann 
weight. Therefore the extraction of the correlations function from the corresponding com- 
mutator Green's function is not as straight forward as in the zero temperature case but the 
overlap to be reversed, i.e. an additional factor cancels out the contribution of the conjugate 



FT((S' A ),(S' B ) 1 h')=0 




(35) 




(36) 
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n 

FIG. 3: Two dimensional density of state 10(71,72) for the bipartite bcc lattice. The invariance of 
the DOS under 71 — > —71 reflects the discrete translational symmetry of the bcc lattice along 
the diagonal, i.e. {x, y, z} — > {x + a/2, y + a/2, 2; + a/2}. The main contribution to the density 
of states comes from a relatively small area around three lines in the shape of the letter Y with a 
peak at 71 = and 72 = —0.327. 

function. This was done by Tyablikov 0, Q in 1959. Here we follow a slightly different 
derivation. 

As we have seen in the disquisition on the zero-temperature regime the adequate self- 
consistency equations are expressed by: 

Wo)*sy> = ~ (37) 

The spectral representation of the Green's function at finite temperature is according to 
Equation (j3J) given by: 



y-v 1 (n\Sf \m)(m\S v j \n) e _/3 ' 



-E 



2n uj — {oJ m — oo n ) He Z 
1 (n\Sj\m)(m\Sf\n) e _/3t " 
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^ 1 (n\Sy \m){m\Sy\n) e _/3a; ™ - e~^ m 
^ 2iru — (uj m — uo n ) ±ie Z 

m,n 

where we once again used that: 

— / e^-^O^dt = lim - 

27T J_ 00 2lX UJ - UJ' + 16 

\m) and \n) are eigenstates of the Hamiltonian and Z refers to the partition function. 
The Fourier transform of the correlation function is given by: 



(38) 



(39) 



<S3W> W = £Hsf|m>Hs»- 

m,n 

x5(uj - (iU m - U) n )) 



-fan 



(40) 



The relation to the Green's function is readily obtained by the following rearrangements: 

m,n l 
e l3(ui n -ui m ) _ j 

\- (n\Sf\n)(n\S]\n) 

n Z 6 

E m , n >(n\Sf\m)(m\S]\n) e -**» - 
e^-l Z 

X 5(U - (uj m - LO n )) 

+ e ~fan 



(41) 



The prime ' on the sum excludes terms where m = n and was introduced to avoid singularities 
in the denominator. The omitted terms are accounted for in the second term of the RHS. 
The following symbolic identity is valid for real uo 

1 1 



lim 

and immediately yields: 



= — 2iii5{u) — uj') 



(42) 



i lim 



ofa _ I 



(43) 
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As pointed out by Stevens and Toombs [13| the second term on the right side is generally 
difficult to calculate. In principle it would have been possible to eliminate this term by also 
calculating the anti-commutator Green's function. Once again, as we have broken the U(l) 
symmetry of the ground state in a way such that (S y ) = this term will yield zero. 
Using the Green's function's Fourier transform into k-space we obtain: 

(Sf(t)S]) = 

i lim / / d 3 k dui 




e^oj " e^-1 

[G yy {k, u + ie)- G yy {k, u - ie)\dw 

fG yy (k,u)e iuJt 



-2tt J d 3 ke ikr ^ [ Residue ^* 



odw _ 1 ' ^ 



+Residue ( — , — uj v 

\ eP» - 1 p / 

(44) 

Now we can define two self-consistency equations which determine the spin fields at finite 
temperature: 

,G^(k,cu) 



F A := J d 3 k ^[Residue(- 



+Residue( - ^ = 



F B := J d 3 k E[ R - d -(^^ } 



+Residue(^^ ) -u, p )] - \ = (45) 



Using the DOS of Equation fl36|) those are written as: 

F A = J djx d^ 2 p{ji, 72) x 

^PD 'A f G A (71,72,^) x 

}jResidue( -^—^ , u p ) 
+Residue( — _ - — , -u p )\ - - = 



F B = I rf7i^72P(7i>72)x 

,Gfl(7i,72,w) 



[Residue (- 



+Residue( G ^ 1 '_ 72 i '^ ) , - Wp )] -1 = (46) 



16 



Again, in the canted anti-ferromagnetic and canted ferromagnetic phases we use the mean- 
field equations to eliminate (S A ) and (S B ) and obtain: 

F%((S%,{SZ),h',T) = 

F c B ((S x A ),(S x B ),h*,T) = (47) 

Similarly, we obtain the self-consistency equations for the ferromagnetic and anti- 
ferromagnetic phases: 

FT((S z A ),(S z B ),h z ,T) = 

F™((S z A ),(S B ),h°,T) = (48) 



V. GREEN'S FUNCTIONS IN EXPLICIT FORM 



In the remaining three Chapters we give the Green's functions and the corresponding cor- 
relations functions and carry out the derivation of the matrix equation defining the Green's 
functions. 

We start with the Green's functions for the anisotropic Heisenberg model for the canted 
phases, i.e. for the canted ferromagnetic and the canted anti-ferromagnetic phases 



-iu(S%) 



■ x 



G ^^ u) (2tt)^-uj 2 1 )(^-uj 2 ) 

{u 2 + ACD + ADA, - 4(C + A 2 )(B 2 + A 2 (A(S X B ) 2 + A(S Z B ) 2 - 1))) 

(lu 2 + ACD + ADA 2 - 4(C + A 1 ){B 1 + A 1 (4(^) 2 + A(S Z A ) 2 - 1))) (49) 



(D(lu 2 + ACD + ADA 2 ) + (B 1 + A 1 (A(S X A ) 2 + A(S Z A ) 2 - 1)) 
(u, 2 - 4(C7 + A 2 ){B 2 + A 2 {A(S B ) 2 + 4(S B ) 2 - 1)))) 

G y B y (k,u) 



(2tt) 4 (^ 2 - toDiu 2 - u 2 2 ) 
{D{u 2 + ACD + ADAi) + (B 1 + A 1 {4,{S X B ) 2 + A(S Z B ) 2 - 1)) 

(u, 2 - 4(C + A^B, + AMSX) 2 + 4(S Z A ) 2 - 1)))) (50) 
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(lu 2 + ACD + ADA, - 4(C + A 2 )(B 2 + A 2 (A(S X B ) 2 + A(S Z B ) 2 - 1))) 

1 • ' ^ j ~ (27r) 4 (^ 2 - w?) (u; 2 - lu 2 ) 
(to 2 + ACD + ADA 2 - 4(C + A 1 )(B 1 + A X {A{S X A ) 2 + A(S Z A ) 2 - 1))) (51) 

The poles are given by: 

u l/2 = 2(A 1 (B 1 + A X (A(S* A ) 2 + A(S Z A } 2 - 1)) + 
A 2 (B 2 + A 2 {A(S B ) 2 + A(S Z B ) 2 - 1))) - ACD 
±2[(-2O0 + A X (B X + A X (A{S A ) 2 + 4(Si> 2 - 1)) 
+A 2 (B 2 + A 2 (A(S X B ) 2 + A(S Z B ) 2 - l))) 2 
-4(C 2 - A X A 2 )(D 2 - (B 1 + A(4(^) 2 + A(S Z A ) 2 - 1)) 
x(fi 2 + A 2 (4(^) 2 + 4(5| ? ) 2 -l)))] 1 / 2 (52) 

All parameter A X} A 2 ,B X ,B 2 , C and D £1X6 cLS stated in Equation (l2Tj) . 

As condition ffTOl) does not apply for the ferromagnetic nor for the anti-ferromagnetic 
states, these Green's functions are structurally different: 

iu{(A 3 + A A )B 3 (S Z B ) - (- 2 oj 2 + A 2 + B 3 B A ){S Z A )) 



G x A y (k,u) 



(2tt) 4 ( 2 u; 2 - 2 ujj)( 2 uj 2 - 2 lv 2 ) 



rxy(h v _ M(^ + A 4 )B 4 (S Z A ) - (- 2 u 2 + A 2 + B 3 B 4 )(S Z B }) 
_ (27r) 4 ( 2 c, 2 - 2 c, 2 )( 2 c, 2 - 2 ^ 2 ) 

G y A y (k,u) ' 



(27r) 4 (cu 2 -^ 2 )(cu 2 -cu 2 ) 
((^ 2 A 3 - A 3 A\ + A 4J B 3J B 4 )(^) + 5 3 (^ 2 + A 3 A 4 - B 3 B 4 )(S Z B )) 



1 



with 



(5 4 (V + A 3 A 4 - J B 3J B 4 )(^1) + (VA 4 - A 2 A 4 + A 3J B 3 S 4 )(5|)) (54) 



A 3 + A 4 ±^(A 3 -A 4 y + AB 3 B A 
ui/2 = o ( 55 ) 
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A, = 2h z + 4(^>(4' - J 2 T 72 (A:)) + ±{S%)Jl 
A, = 2h z + A(S Z B )(4 - J 2 T 7 2 (A:)) + A(S Z A )4 

B 3 = -A(SX)Jjli(k) 

B 4 = -MS^jJ^ik) (56) 



A. The Correlation Functions in Explicit Form 

Here we state the relevant correlation functions. All other correlation functions yield 
identically \{S XA/B ) or \{S ZA/B ) and are therefore not suitable as self-consistency equations. 
We start with the Green's function for the canted anti-ferromagnetic and the canted ferro- 
magnetic phase at finite temperature: 



(aCD 2 + u)\D + u\B x - ACB X B 2 
-AAlB^Sl) 2 + A(S Z " 2 1A 



Bl 

+A 1 (A(S X A ) 2 + A(S Z A ) 2 -1) x 
(col - 4CB 2 - AA 2 2 {4(S X B ) 2 + A(S Z B ) 2 - 1) 
-AA 2 (B 2 + C(A(S B ) 2 + A(S Z B ) 2 - 1))) 
+AA 2 (D 2 - B X {B 2 + C(A(S X B ) 2 + A(S Z B ) 2 - 1)))) 

coth(^) / 
- of 3 2 \ ( ACD + "l D + ufa - ACB X B 2 

-AAjB^S^) 2 + A(S B ) 2 - 1) 

+A 1 (A(S X A ) 2 + A(S Z A ) 2 -1) x 

(uj 2 2 - ACB 2 - 4A 2 (4(^> 2 + 4(^> 2 - 1) 

-AA 2 (B 2 + C(4(^) 2 + 4(^> 2 - 1))) 

+AA 2 (D 2 - B,{B 2 + C(A(S X B } 2 + A(S Z B } 2 - 1)))) (57) 

(S*s> (Q)), cothm x 

(4C/J 2 + cj 2 /} + c^£ 2 - 4C525! 
-4A 2 5 2 (4(^) 2 + 4(^) 2 -l) 
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+A 2 (4(^) 2 + 4(^> 2 -l)x 
H 2 - ACB 1 - 4A 2 (4(^) 2 + 4(S Z A ) 2 - 1) 
-AA^B, + C(4(^) 2 + 4(^) 2 - 1))) + 
AA^D 2 - B 2 (B 1 + C(A(S X A } 2 + A(S Z A } 2 - 1)))) 

- 7-7 3 V 2 \ (ACD + u 2 2 D + lu 2 B 2 - ACB 2 B X 

-4A 2 J B 2 (4(^) 2 + 4(^) 2 -l) 

+A 1 (4(^) 2 + 4(^) 2 -l)x 

(lo 2 2 - ACB 1 - AA\(A(S X A ) 2 + A{S Z A } 2 - 1) 

-AA^B, + C(4(^) 2 + 4(^) 2 - 1))) 

+AA 1 (D 2 - B 2 (B 1 + C(A(S X A ) 2 + A(S Z A ) 2 - 1)))) (58) 



At zero temperature these functions become: 

W^(0))> = - o * ,„ x 



(4CT> 2 + cj 2 L> + c^Bi - ACB X B 2 

-4A 2 5 1 (4(^) 2 + 4(^) 2 -l) 

+A 1 (4(^) 2 + 4(^) 2 -l) x 

K 2 - 4C5 2 - 4A 2 (4(S|> 2 + A{S%) 2 - 1) 

-AA 2 {B 2 + C(A(S X B ) 2 + A(S Z B ) 2 - 1))) 

+AA 2 (D 2 - B X {B 2 + C(A(S X B ) 2 + A(S Z B ) 2 - 1)))) 

UCD 2 + u 2 2 D + cj 2 ^ - ACB X B 2 



V W 2( W 2- W l) V 

-4A 2 5 1 (4(^) 2 + 4(5| ; ) 2 -1) 
+A 1 (4(^) 2 + 4(^) 2 -l) x 
(u 2 2 - ACB 2 - AA 2 2 (A(S X ) 2 + 4(^> 2 - 1) 
-AA 2 {B 2 + C(A{S X B ) 2 + A(S Z B ) 2 - 1))) 

+4A 2 (D 2 - B 1 (B 2 + C(A(S X B ) 2 + 4(5'| J > 2 - 1)))) (59) 



2^K 2 -a; 2 ) 
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(aCD 2 + uj\D + u\B 2 - ACB 2 B 1 
-AA\B 2 (A(S X A ) 2 + A(S Z A ) 2 - 1) 
+A 2 (4(^) 2 + 4(^) 2 -l)x 
{ul - ACB 1 - AAlmS'X) 2 + A(S Z A ) 2 - 1) 
-AA^B, + C(A(S X A } 2 + A(S Z A } 2 - 1))) + 
AA 1 (D 2 - B 2 (B X + C(A(S X A ) 2 + A(S Z A ) 2 - 1)))) 

=1 (4O0 2 + cj 2 L> + w|S 2 - 

2v^(c 2 2 -c 2 )V 

-4A 2 5 2 (4(^) 2 + 4(^) 2 -l) 

+A 1 (A(S x B ) 2 + A(S z B ) 2 -l) x 

(uj 2 2 - ACBi - AAl(A(S x A } 2 + A(S Z A } 2 - 1) 

-AA 1 (B 1 + C(A(S X A ) 2 + 4( 1 S1> 2 - 1))) 

+AA 1 (D 2 - B 2 (B, + C(A(S X A ) 2 + 4(^> 2 - 1)))) 

the correlation functions for the ferromagnetic and the anti-ferromagnetic phases 

(SfS] eA (0)) k = 

coth(^) x 
2(uf - uxuf) 

((A 3 (S A ) + B,{S z B ))ul - (A 3 A 3 - B 3 B,)(A 4 (S A ) - B 3 (S Z B ))) 
, coth(^) x 
2(u;f - u; 2 cu 2 ) 

((A 3 (S A ) + 5 3 (^»^ 2 2 - (A 3 A 4 - B 3 B A )(A A (S Z A ) - B 3 (S Z B })) 

(SfS] eB (0)) k = 
coth(^f) 
2(uj\ - ujxujIY X 

((A A {S Z B ) + Bi{S z A ))Jt - (A A A A - B A B 3 )(A 3 (S Z B ) - B 4 (S Z A ))) 
coth(^) 

_| __±_1 x 

2(a>2 — ujfu 2 ) 

((MS'b) + 5 4 (^))c 2 - (A 4 A 3 - B 4J B 3 )(A 3 (^) - B A (S Z A })) 
which can be written in the zero temperature limit as: 



(s?sy(v)h = 
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(A 3 A A - B 3 B 4 )(A 3 (S ZB ) - B A (S ZA )) + (A 4 (S ZB ) - B 4 (S ZA )) y^^F (g4) 



B. Derivation of the Matrix Equation 

Finally we carry out the derivation of the matrix equation, evolving from the equations 
of motion. The three equations of motion for the Green's functions after employing the 
cumulant decoupling scheme become: 

id t (%(t)=i5(t)8v(sm)-ih,<%(t) 

-2iJ24((sf(t))G z y(t) + (smog®) 

+2* Ju«Sm)<%(t) + (S?(t))G*](t)) (65) 
i 

id t G\]{t) = ih z G^{t) 

+2iJ24((Sm)G%(t) + (SKt))G%(t)) 
I 

-2. Jl((S-(t))Gl y (t) + {S?(t))G%(t)) (66) 
i 

id t G%(t) = -i5(t)8ij(S?(t)) 

+2iJ2Ju((S-!(t))G%(t) + (Sf(t))G*(t)) 
i 

-2i JZ((Sf(t))G%(t) + (S?(t))G%(t)) (67) 
i 

Fourier transforming into cj-space and assume that the spins are constants in time: 

u(%(u) = ^6 ij {Sl)-ih z C%(u) 

-2iY,4{(WG%(u>) + (S;)G!g(u>)) 
i 

+2* J7i((S?)G™(u,) + (Sf)GtJ(co)) (68) 
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uCg(u) = ih z G*]{u) 
i 

-2. Jj((S?)C%(u) + (S?)G%(u)) (69) 
i 

+2 ? ^jJ((5f)G^H + (^)G^M) 

-2. £ J^((S*)G™(u) + (SftGgV)) (70) 

In the canted anti-ferromagnetic and the anti-ferromagnetic phases the spins on the two 
sub-lattices assume different values. Therefore we need to split up every Green's function 
before we can Fourier Transform into k-space. We define G A ^(k, uj) and G B u (k, uj). G^i^k, uj) 
refers to G^ {uj) when the site i is on sub-lattice A and G B u (k, uj) when i is on sub-lattice B. 

-4i(S? A )[G A ?(u)jh2(k) + G B ?(u)jhi(k)] 
-4<G^(a;)[(5i>4 l 72(0) + (S' B )Jhi(0)] 
+MS z A }[G A f(uj)jJ l2 (k) + Gsf^jJ^k)] 
+4iG A ? (uj)[(S y A )jJ l2 (0) + (S B ) 7771(0)] 

(71) 



uGjg(u) = -2ih,G A ?(u) 

+4i(S° A )[G A ?(u)4>y2(k) + G B ?(u)jh(k)] 
+4iG A ? (u)[{S z A ) J| 72 (0) + (S' B ) Jhi(0)} 
-Ai(S\)[G A ^{u)J^ l2 {k) + G B ? (uj) jj^(k)\ 
-UG A ? (u)[(S* A ) J 2 T 72 (0) + (^)J 1 T 7 i(0)] 

(72) 

uG A ? (uj) = -j^iiSZ) 

+MS v A )[G A x k y (uj)Jj l2 (k) + G B ?(u>)J?>n(k)] 
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-AiG A f {u)[(SD J 2 T 72 (0) + (^)J 1 T 7 i(0)] 
-±i(Sl)[G A f (u) Jj l2 (k) + Gsf^jJ^ik)] 
-4*G4>)[(^)J 2 T 72 (0) + (SI) 7771(0)] (73) 



VI. CONCLUSIONS 



In conclusion, we analysed the anisotropic Heisenberg model in a external field on 
the three dimensional bcc lattice by employing the well-established technique of real-time 
Green's functions for spin systems. The series of infinite order Green's functions as it ap- 
pears in the equation of motion was truncated by applying cumulant decoupling and the 
resulting random phase approximation accounts for linear spin-waves. We are the first to 
apply this method to the canted anti-ferromagnetic phase entailing a set of six algebraic 
equations. The innate self-consistency equations herein constitute a three dimensional nu- 
merical integral over the k-space. By introducing a two dimensional density of states the 
integral was reduced to two dimensions where the lattice generating functions serve as new 
integration variables. In the appearing integrals the DOS is the only quantity that depends 
on the structure of the lattice. Hence, once the DOS is computed for a certain lattice ge- 
ometry the further calculation remain unaltered. Therefore our method is widely applicable 
and easily adjustable to various magnetic systems where canted phases are in the center of 
interest. This also holds for two dimensional lattices where linear spin waves are expected 
to yield a reasonable approximation. 
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APPENDIX A: CONNECTION TO 4 HE PHYSICS 



For the interested reader, we briefly summarize the connection between the anisotropic 

nn 

Heisenberg model and He as introduced by Matsubara and Matsuda [l4j, [15| , and used by 
Matsuda and Tsuneto [16| . Fisher 5J, Liu and Fisher |y] and most recently by Stoffel and 
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Gulacsi 0,13,3. 

Apart from possible 3 He impurities 4 He is a bosonic system and the generic Hamiltonian 
for such systems in the language of second quantization is given by: 



H = I d 3 xfthc)( — -V 2 + /i)V>(x) 
J 2m 



(Al) 



where ^(x), the particle creation operator and ^(x) , the corresponding destruction oper- 
ator obey the usual bosonic commutator relations. Hamiltonians in three dimensions such 
as in Equation flAlj) are not solvable even for elementary potentials V(x) such as the Dirac 
delta distribution. Therefore, further approximations has to be implemented. An approxi- 
mation which proved particularly successful for the description of liquid Helium is know as 
the quantum lattice gas model and was first introduced by Matsubara and Matsuda 

In the quantum lattice gas model one works with a space lattice of discrete lattice points 
rather than the continuum. This approximation proves to be very useful to study solid states 
as the spacial discretization of this model serves as a natural frame for the crystal lattice. 
Also in this model no specific knowledge of the density distribution of the atoms is needed. 

According to Matsubara and Tsuneto [16] the generic Hamiltonian Equation (1A1I) in the 
discrete lattice model reads: 

H = n 2J rii + 2J Uij{a\ - a]){di - a,j) + ^ VijUiUj 



(A2) 



Here Uij are non-zero for nearest neighbor and next nearest neighbor hopping and otherwise 
zero. The values of u nn and u nnn are such that the kinetic energy is isotropic up to the 4th 
order. In the case of a bcc lattice (two interpenetrating sc lattices) the matrix elements are 
given by: 

Unn = I w (A3) 

Unnn = -~. 3 ( A4 ) 

As the atoms do not penetrate each other there can exist only one atom at a time on a 
lattice site. Consequently and a are the creation and annihilation operators of a hard 
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core boson commuting on different lattice sites: 



[a>i, «jj 



but obey the anti-commutator relations on identical sites: 



(A5) 




1 



(A6) 



Equation (1A2jl is the Bose-Hubbard model in three dimensions for hard core bosons. Due 
to the unusual statistics of hard core bosons, Wick theorem cannot be applied and hence, 
the common formalism of perturbative field theory is not applicable. The way out is to 



transform the model to an equivalent spin model [14 

at = Sf - is* 



151 ]. namely by using 



71; 



S7 + z5- 
- 2 *j 



(A7) 



It can be verified that the usual Lie algebra for spin 1/2 particles preserves the mixed 
commutation/anti-commutation relations for hard-core bosons. This substitution trans- 
forms the hard-core bosonic Hubbard model into a spin model: 

i 

+ E - s " - s ] - 2S " S ! - 2S ! s j) 



H 



+E^4-f -f +s & 



(A8) 



If we adjust the notation to conform with the usual standards of spin models, we re-obtain 
the anisotropic Heisenberg from Equation ([T]): 



(A9) 



ij 



with: 



J 



J 



—2ui 



-» + E J l - E 4 



(A10) 
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If the above presented transformation is used for 4 He, then the values of the J's also have 
to be chosen such as to mimic 4 He. The interactions between the 4 He atoms are controlled by 
van-der-Waals forces and their repulsive nature at very short distances determines negative 
nearest neighbor interaction jf, evoking anti- ferromagnetic ordering in the spin language. 
The corresponding Lennard- Jones potential is short ranged and therefore it is sufficient 6| 
to only consider nearest and next nearest neighbor interactions. Hence, for 4 He the J values 

will be if = -qij} eAjeB , 4 = ^jfeAjeAi J 7 = ~Qi J 7eAjeB and J J = -^JJzAjeA, where 
qi = 6 and q<i = 8 are the number of nearest and next nearest neighbors on the bipartite bcc 
lattice. Liu and Fisher 6] in their calculations used, J7 = 1AK, Jj = 0.5K, jl = — 3.8-K" 

Tin 

and = —1.7K. In a 4 He calculation the results do not change |3|, |4J if the J"'s are within 
±2 range of these values and J B > jjland J T 's v&. 



Defining two sub-lattices gave 



14 
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ues remain positive. 



16| a possibility to establish the diagonal 
long-range order of solids in a natural way: sub-lattice A represents the centers of the 4 He 
ions, hence it coincides with the ion lattice. Sub-lattice B defines the interstitials, the space 
in-between those atomic centers. In the liquid phases, of course, the occupation number 
on both sub-lattices is equal as there is no spacial density variation. In Table (111) we gave 
the various magnetic phases of the anisotropic Heisenberg model. These phases, however 
identify the corresponding phases of the 4 He system, as presented in Table ( llT|) . 



APPENDIX B: THE MEAN-FIELD LIMIT 



In the mean-field solution of Liu and Fisher [6( the Green's functions have not been 
explicitly evaluated, for the sake of completeness we re-derive these hereafter. 

First we will re-derive the classical mean-field approximation as was pioneered by Liu and 
Fisher js] and briefly state some key properties. We will further show that this approximation 
is a special case of the random-phase approximation. 

The anisotropic Heisenberg Hamiltonian in the classical mean-field approximation is ob- 
tained by substituting the spin 1/2 operators with their respective expectation values: 

H M f = -h z {{S\) + (S' B )) 

-24(SX)(S B ) - 4((S Z A )(S Z A ) + (S%)(S%)) 
-2J?(S%(S%) - J 2 T «Si)<^> + (S%)(S%)) 
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Spin Configuration 


Magnetic Phase 


ODLRO 


DLRO 


4 He-Phase 


TT 


FE 


No 


No 


Normal Liquid 


// 


CFE 


Yes 


No 


Superfluid 


// 


CAF 


Yes 


Yes 


Supersolid 


U 


AF 


No 


Yes 


Normal Solid 



TABLE II: The phases of A He corresponding to the phases of the anisotropic Heisenberg model. 
Similarly to Table ((J) the phases are defined by their long range order, i.e., off-diagonal long-range 
order (ODLRO) and diagonal long-range order (DLRO). 

(Bl) 

Here j\ = -qij! eAJeB , 4 = -92./jeAje^ J 7 = -qiJjeAjeB and J J = -^JJeA^A where 
qi = 6 and q^ = 8 are the number of nearest and next nearest neighbours on the bipartite 
bcc lattice. The mean value of S y drops out as the randomly broken symmetry S x <-» S y 
(off-diagonal long-range order) allows for (S y ) = 0. The standard method of deriving the 
corresponding self-consistency equations is to minimize the Helmholtz's Free energy F = 
H — TS. The entropy S is given by the pseudo spin entropy of the system: 

S = ~\[{\ + S A ) ln(il + S A ) + (i - S A ) \n(\ ~ Sa) 

+(i + S B ) ln(i + S B ) + (5 - S B ) In(i - S B )] (B2) 

where Sa = a/ (S za } 2 + (S XA ) 2 and S B = \J (S ZB ) 2 + (S XB ) 2 . We could equally well say that 
the state of the system is determined by minimizing the internal energy (H), subject to an 
additional constraint given by Equation (1B2I) . In this picture the temperature becomes a 
Lagrange multiplier and at absolute zero, where T = S = we obtain, as lim5 A s B ^ S = 0, 

\j(siy + (si) 2 = ^{s B ) 2 + {s B Y= l - (B3) 
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This result implies that this approximation does not take quantum fluctuations into account. 

In the canted anti-ferromagnetic and the canted ferromagnetic states there are four self- 
consistency equations in the ferromagnetic and anti-ferromagnetic phases; where (S x ) = 
(S y ) = they are reduced in number by two. These equations are readily obtained by 
differentiating the free energy with respect to (S z ) and (S x ) respectively. The resulting 
equations can be rearranged to yield: 

2tu A 

(S Z A ) = 2J l'(^) + 2jJ(^) + ^ tanh( ^ ) 
2uj a 

/qx , 2Jj(S A ) + 2Jj(S%) 

(S B ) = tanh(/3w s ) 

ZUJ B 

(S B ) = 2jl h S A) + ^4(S z B ) + h z tanH(3ujB) (B4) 
2uj b 

where 

u; A =[(2J 1 T (^)+2J 2 T (^)) 2 + 

{2jl(S B ) + 2jl(S\) + h z f^ 
co B = [(2jJ(S A }+2jJ(S B }) 2 + 

{2j\\(S B )+2j\\(S z A ) + h z )^ (B5) 

In the canted phases where the transversal magnetic fields (S A ) and (S B ) are non-zero, the 
energies u A and uj b can be eliminated from equations (IB4I) to yield the following important 
relations: 

h* + 2(S A )(4 ~ Jj) + 2(S Z B )4 = 2JJ { ^{S A ) 
h z + 2(S Z B )(4 - Jj) + 2(S z A )jf = 2jj r M<S£) 



(B6) 



In the limit h z — > oo the Hamiltonian (Equation (IBlj) ) reduces to an effective single 
operator model: 

H = -h*{(S A ) + (S%)) (B7) 

Consequently for high external fields the system will assume the energetically favorable fer- 
romagnetic phase. In the opposite limit h z —>■ and with sufficiently large anti-ferromagnetic 
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nearest neighbour coupling J{ <C the system is dominated by: 

H = -4(S Z A )(S Z B ) (B8) 

giving rise to the anti-ferromagnetic state. At medium large fields h z the two terms Equation 
(]B7j) and Equation (IB8I) balance each other and the transversal ferromagnetic coupling 
( Jj > and Jj > 0) become significant, deviating the spins into the x-direction. In regions 
of higher h z , where ferromagnetism (Equation ( 1B7I) ) is more prevalent the system leaps into 
the canted ferromagnetic phase; for lower h z it yields the canted anti-ferromagnetic phase. 
For those sets of coupling constant where all four phases are existent, the corresponding 
phase transitions are of second order. If, due to choice of constants one or more of those 
phases, for example the canted anti-ferromagnetic phase, does not appear the resulting 
canted ferromagnetic to ferromagnetic phase transition is first order. 

Now, we may ask the question "For which sets of parameters all four phases appear?". 
Matsuda and Tsuneto Q derived Prices for aa phase rian S1 rion points at zero te m pe r - 
ature. As mentioned in the previously the four phases are distinguished by their order 
parameters, m\ = (S^) + (S B ) and m 2 = (S%) — (S B ). Across all second order phase 
transitions the spin mean-fields (S x ) and (S z ) are continuous C° functions of h z . Therefore 
the critical points are determined by Equations ( 1B6I) . in the limits where the relevant order 
parameter, mi or 777.2, disappears. The canted ferromagnetic phase transits into the ferro- 
magnetic phase when mi approaches zero. Hence we set (S^) = (S B ) ^0 and consequently 
(Sa) = (S B ) = \. Equation (1B6I) readily gives: 

^fe-cfe = J\ + Jj ~ J\ ~ ^2 (B9) 

Equally the canted anti-ferromagnetic to anti-ferromagnetic transition is defined by (S^) = 
— (S B ) = \ while (S^) — > and (S B ) — > 0. The unknown quotient ||#y is readily eliminated: 

h Z C AF-AF = \j (-4 + 4 - JJ) 2 ~ {JJ? (BIO) 

The canted ferromagnetic and the canted anti-ferromagnetic phases coexist where the or- 
der parameter of the diagonal long-range order, 777.2 = (5*1) — (S B ) approaches zero. We 
replace (S A ) and (S B ) in equation (1B6j) with m 2 and m\ and retain only linear terms of m 2 . 
Subtracting and summing up both equations respectively yields: 

h z + 2m2 (4 1 - Jj + j\) = 1Jjm 2 
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2m 1 (4 l - Jj - jf) = -2jJml Am l + l (Bll) 
v 2 2 u 1 Ami - 1 



We used that a/ (S^) 2 + (^l) 2 = 4 at T=0. The solution of these two equations determine 
the critical point which is given by: 

7,2 _ ^1 + J 2 ~ J I ~ Jj 

I^CFE-CAF ~ || _ || - >< 
t/j ^2 J l t o<2 



'(-Jf + Ji 1 -^) 2 -^) 2 (B12) 
For a particular choice of coupling constants all four phases will exists when: 

h Z FE-CFE > hcFE-CAF > ^CAF-AF (B13) 

In other cases, for example where h z CFE _ CAF < h z CAF _ AF the canted anti-ferromagnetic 
phase is suppressed. The resulting first order canted ferromagnetic - anti-ferromagnetic 
transition point has to be calculated by finding the state with the lowest internal energy 
(H) and making a comparison between the two. 



APPENDIX C: LINK TO MEAN-FIELD SOLUTION 

Here we establish a link between the classical mean-field approximation and the random- 
phase approximation as derived in the previous Chapters. We have already seen that both 
approximations are mean-field type, involving self consistency equations and that there are 
two equations (Equation flB6j) or Equation (j!9p ) which appear in both approximations. 
However, the random-phase approximation takes spin-wave/quasi-particle excitations into 
account whereas the classical mean-field approximation is an effective one operator model 
exhibiting two energy levels per sub-lattice. Therefore we can now consider if those en- 
ergy levels, specifically their difference, correspond to certain spin-wave excitations. In 
this Appendix will show that the classical mean-field approximation is a limiting case of the 
random-phase approximation in a way in which the integral of the momentum k is restricted 
to values where: 

7i(*0 = 72(*) = 0, (CI) 
or equally the generalized density of states is: 

p(7i,72) = 5(7i)% 2 ) (C2) 
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FIG. 4: The rings show the areas where the two lattice generating functions 71 (A;) and 72 (k) yield 
simultaneously zero. If the integral over the first Brillouin zone is restricted to those areas the 
random-phase approximation devolves into the classical mean-field approximation. 



In this limit the self-consistency equation of the random-phase approximation Equation ( 146]) 
becomes 



f a M f = [Residue ( 



Gff(7i = 0,72 = 0, op 



, G^(7l = 0,72 = 0,a;) » 1 n 

+Residue( _ , -u p )\ - - = 

f ^rp ^ , ^(71 = 0,72 = 0,^) , 

^ , Gf(7i = 0,72 = 0,^) v, 1 n ™ 
+Residue( _ , -u p )\ - - = (C3) 

Figure H] shows the wave- vectors that correspond to 'ji(k) = 72 (A;) = within the first 
Brillouin zone. Physically, vanishing lattice generating functions (71 (A;) = ^(k ) = 0) means 
that the system loses all information about the lattice structure as it is the case in the 
classical mean-field approximation where the only information that remains is the number 
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of nearest and next nearest neighbors. With the lattice generation functions being zero the 
matrix M of Equation ffl2l) becomes: 

/ 








Mi 3 














%lo 





M 24 








M 13 





iuj 





M 35 








-M 24 





za; 





M 46 








M 53 





iu; 














M 64 








M 1 



V 

where the components are given by: 

M 13 = 2h z + 4(S Z A )4+4(S Z B )J 1 } 
M 24 = 2h z + 4(S Z B )4 + H$a) 4 
M 35 = A(S X A )JJ + A(S B )JJ 
M 46 = -4(^)J 2 T +MS%)JJ 
M 53 = -4(S%)JJ-4(SI)JJ 
M 64 = -4(SI)JJ-A(S X B )JJ 

The corresponding Green's functions are given by: 

(AJJ(S B )+UJ(S%))(S X A ) + (44(S z B )+44(S z A ) + 2h z )(S z A ) 



(C4) 



(C5) 



GT (w) 



CO 



(2^0 



J) 



(AJ^(S x A )+AJj(S xB ))(S' B ) + (AJl{S z A ) + A4{S z B ) + 2h z ){S z B ) 



u 2 - (2uj 2 y 

The poles of the Green's functions are given by the eigenvalues of M cmf : 

u x = ±[{2Jj(S%) + 2JJ (S x A )u A ) 2 + 
(24(S z B )+24(S z A ) + h z )^ 

uj 2 = ±[{2Jj(S x A ) + 2J 2 T (S B )u A ) 2 + 
(2j\S z B ) + 2jW{S z A ) + h z ) 2 }^ 



(C6) 



(C7) 



Those energies are, as we have expected, identical to the classical mean field energies given 
by Equation (1B5I) . The self-consistency equations calculated with Equation ( 1C3I) yield: 

1 



2wi 



[{2Jl (S%) + 2J*(S* A )){S* A ) + 
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{2jl(S* B )+24(S'jd + h')(S' A )] 
1 



tanh(/3u;) 



2v< 



\{2Jl (S A ) + 2J 2 (S xB ))(S* B ) + 



(24(S A )+24(S B ) + h*)(S B )} 



tanh(/?a;) 



(C8) 



(C9) 



Those two equations together with Equation (TT9l) are identically to Equation (1B4I) . Therefore 
we have shown that the mean-field as derived by Fisher and Liu [(3] is indeed a limiting case 
of the random-phase approximation. 
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